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ABSTRACT 

In several GRBs afterglows, rapid temporal decay is observed which is inconsistent with spherical (isotropic) 
blast-wave models. In particular, GRB 980519 had the most rapidly fading of the well-documented GRB after- 
glows, with f- 2 05 ± 04 in optical as well as in X-rays. We show that such temporal decay is more consistent with 
the evolution of a jet after it slows down and spreads laterally, for which t~ p decay is expected (where p is the 
index of the electron energy distribution). Such a beaming model would relax the energy requirements on some 
of the more extreme GRBs by a factor of several hundreds. It is likely that a large fraction of the weak (or no) af- 
terglow observations are also due to the common occurrence of beaming in GRBs, and that their jets have already 
transitioned to the spreading phase before the first afterglow observations were made. With this interpretation, a 
universal value of p = 2.5 is consistent with all data. 
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1. INTRODUCTION 

One of the most important open questions in GRBs is 
whether the burst emission is isotropic or strongly beamed in 
our direction. This question has implications on almost every 
aspect of the phenomenon, from the energetics of the events, 
to the engineering of the "inner engine" to the statistics and 
the luminosity function of the sources. We suggest here that 
GRB 980519 was a jet with an opening angle of less than 0.1 
rad. We also suggests that such jets are common in most GRBs. 

According to the relativistic fireball model, the emission 
from a spherically expanding shell and a jet would be rather 
similar to each other as long as we are along the jet's axis and 
the Lorentz factor, 7, is large compared to the inverse of the 
angular width of the jet, 60 (Piran, 1995). When 7 drops below 
0Q 1 the jet's material begin to spread sideways and we expect a 
break in the light curve of the afterglow at this stage. Since we 
have for spherical adiabatic evolution 7(f) w 6(Es2/n\) l ^t d ^ 8 , 
this break should take place at 3 



knowledge of the value of the electrons' energy distribution in- 
dex p ~ 2.4 we expect for high frequencies a spherical decay of 
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where £52 is the "isotropic" energy of the ejecta in units of 
10 52 ergs, i.e., the inferred energy assuming isotropic expansion, 
and n\ is the surrounding ISM particle density in cm" 3 . So far, 
with the exception of the recent GRB 990123 (Kulkarni et al. 
1999), no such break was observed, even for afterglows extend- 
ing for hundred days. More specifically, the best observed after- 
glows, those of GRB 970228 and GRB 970508, behave accord- 
ing to a single unbroken power law, as long as the observations 
continued (Zharikov, Sokolov, & Baryshev 1998; Fruchter et 
al. 1998), giving a strong indications that those sources were 
isotropic to a large extent. 

We show here that, even without seeing a break in the 
lightcurve, one can identify a jet based on the powerlaw index 
of the optical light curve decline. Since we have a reasonable 
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and a jet like decay of t . We suggest that at least in 
one afterglow, GRB 980519, the observed light curve and spec- 
tra are consistent with an expanding jet and inconsistent with 
those expected from a spherical expansion. We suggest that in 
this burst the transition to spreading jet, at 7 ~ Oq 1 , took place 
during the few hours between the GRB observations and the 
first detection of the afterglow. We conclude that the beaming 
factor in this burst is at least a few hundred. Together with the 
appearance of a sharp break in the light curve of the afterglow of 
GRB 990123, this indicates that jets are common in GRBs. In 
fact the rapid decline that corresponds to an expanding jet could 
also explain the weak or no optical afterglow seen in some of 
the other bursts, e.g., GRB 990217 (Piro et al. 1999; Palazzi et 
al. 1999). 

Jets have been discussed extensively in the context of GRBs. 
First the similarity between some of the observed features of 
blazars and AGNs led to the speculation that jets appears also in 
GRBs (Paczynski 1993; Dermer & Chiang 1998). Second, the 
regions emitting the GRB as well as the afterglow must be mov- 
ing relativistically. The emitted radiation is strongly beamed 
and we can observe only a region with an opening angle 1 /j 
off the line of sight. Emission outside of this very narrow cone 
is not observed. This have lead to numerous speculations on the 
existence of jets and to attempts to search for observational sig- 
nature of jets both during the GRB phase (Mao and Yi 1994) 
and in the context of the afterglow (Rhoads, 1997a,b;1998; 
Meszaros et al., 1998; Panaitescu & Meszaros 1998). Finally, 
GRBs appear naturally in the context of several leading scenar- 
ios for the "inner engine". (Mochkovich et al., 1993; Davies 
et al., 1994; Katz, 1997, Meszaros & Rees, 1997, Nakamura, 
1997). 
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and Panaitescu, & Meszfos (1998). We explain these differences in section 2. 
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The simple fireball model [and the Blandford-McKee (1977) 
solution] assumes a spherical expansion. However, even if the 
relativistic ejecta is beamed then as long as the Lorentz factor 
7 of the relativistic motion satisfies 7 > , the hydrodynam- 
ics of the jet won't be influenced by the fact that it has a finite 
angular size (Piran, 1995). The matter doesn't have enough 
time (in its own rest frame) to expand sideways. This situation 
changes drastically when 7 w Oq 1 when the sideway expansion 
becomes significant. A full solution of the evolution of a jet at 
this stage requires 2D relativistic hydro simulations. However, 
one can obtain a reasonable idea on what goes on using simple 
analytical estimates. 

Pvhoads (1997a,b;1998) considered the evolution of a rela- 
tivistic jet that is expanding sideways at the local speed of 
sound, c.5 so that 9 ~ 6$ + c s t proper / ct ~ 9q + 7" 1 / v3. In this 
case the hydrodynamic transition takes place at 7 ~ 9~ l /V3. 
However, as the rest mass of the shocked material is negligi- 
ble compared with its internal energy, the expansion can be ul- 
tra relativistic with a Lorentz factor comparable to the thermal 
Lorentz factor. This would lead to 9 ~ 9o+ct proper /ct ~ 9q + j~ 1 
and to a transition when 7 ~ 9q . The sideways expansion leads, 
for an adiabatic evolution, to an exponential slowing down as 
7 oc exp [—r/ljet], where lj e , = [Ej et /(4Tr/3)nm p c 2 ] 1 ^ is the Se- 
dov length in which a spherical expanding shell with energy 
Ej et acquire mass whose rest mass energy equals to its own en- 
ergy (n is the ISM density). Ej el is the actual energy in the 
jet. Thus, r is practically a constant during the spreading phase. 
Therefore, the observer time, which is related to the radius and 
the Lorentz factor as t oc r/7 2 satisfies simply t oc 7 ~ 2 . 

Our estimate for the break time (Equation [l]) is the simplest 
one. It is based just on spherical adiabatic expansion. It dif- 
fers by a factor of 20 in time (corresponding to a factor of 
~ 3 in the opening angle #0) from the expression given by 
Rhoads (1998). The discrepancy arises from several factors: 
(i) As discussed above, we assume that the jet expands side- 
ways at the speed of light while Pvhoads (1998) assumes that 
jet expands at the sound speed c/v3. (ii) Rhoads (1998) uses 
t = R/2j 2 c. This expression is valid for a point source mov- 
ing along the line of sight with a constant velocity. We use 
t w R/4^f 2 c reflecting the deceleration of the source and its fi- 
nite angular size (Sari, 1997,1998; Waxman 1997, Panaitescu 
and Meszaios 1998). (iii) We use the simple adiabatic energy 
condition: E = ^ 2 mc 2 , where m is the rest mass of the shocked 
ISM, while Rhoads (1998) uses E = 2 7 2 mc 2 . A third possibil- 
ity is to use the more exact numerical factor derived from the 
Blandford McKee (1976) solution: E = 12 7 2 mc 2 /17. (iv) We 
estimated the time in the local frame as R/jc. Rhoads noted 
that the Lorentz factor was higher earlier and hence the effec- 
tive proper time is shorter by a factor of 2.5 allowing for less 
spreading. However, far from the shock, the matter moves with 
a considerably lower Lorentz factor allowing it to spread more 
easily. 

Panaitescu & Meszaros (1998) consider similar hydrody- 
namics as Rhoads (1998) but notice that once 7 ~ 1/0q the 
observer is able to see the edge of the jet. They find two tran- 
sitions, the first one when 7 ~ 1 /9q at around our break time 
estimate and the second one around Rhoads'. However, there 
would be only one transition if the time between the two breaks 
turns out to be very short. A reliable estimate of the numeri- 
cal factor clearly requires full 2D simulations. It might also be, 
as suggested by Rhoads( 1998b), that the transition takes place 
over a relatively long time and that most observations, that are 



conducted in a finite time interval, will show only part of the 
asymptotic break. 

We consider now synchrotron emission from a powerlaw dis- 
tribution of accelerated electrons produced by shocks in an ex- 
panding jet. The instantaneous spectrum is given by the four 
broken power laws discussed in Sari, Piran & Narayan (1998). 
However, the time dependence of the break frequency and the 
overall normalization depend strongly on the hydrodynamic 
evolution. Therefore, the lightcurve from a jet differs strongly 
from the light curve of a spherical evolution. Surprisingly, it 
is possible to obtain general expressions, appropriate to both 
spherical and jet evolution (by spherical we mean any system 
with 7 > 9~ l and by a jet a system with 7 < 9~ l ). We write 
these generalized expressions and specialize to jet and sphere 
only at the very end. We begin with the typical frequency v m , 
at the observer frame: 
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The cooling frequency, the synchrotron frequency of electrons 
that cool on the dynamical time of the system, is given by 
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The peak flux is obtained at the lowest of the two frequencies 
v m and v c . Let N e be the total number of electrons radiating 
towards the observer, i.e., those located in a cone of opening 
angle 7" 1 . [N e is different from N e (Sari et al., 1998; Sari & 
Piran, 1999) which is the total number of radiating electrons, 
including those that are not radiating towards the observer]. N e 
can be approximated by N e = ir~f~ 2 R 3 n/3. The total energy 
per unit time per unit frequency emitted by these electrons, 
ajm e c 2 N e B^i /6ne, is distributed over an area of irj~ 2 d 2 at a 
distance d from the source. The observed peak flux density is 
therefore 
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It seems to hold quite generally at late times (except per- 
haps the first few hours, see Sari & Piran, 1999) that v c ^> v m . 
The electrons responsible for low energy emission are therefore 
those with v m . In this case, the self absorption frequency can be 
estimated as 
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We now turn to calculate the light curves for several fre- 
quency ranges. The flux at low frequencies, which is self ab- 
sorbed, evolves as 



F v<v ocfl 2 oc<j fl/2 spherical, 
1 const, jet. 



(6) 



The flux at frequencies that are above the self absorption fre- 
quency but below the typical frequency v m evolve as 
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We therefore expect that the low frequencies (v < v m ) flux 
would rise like t 1//2 as long as the evolution is spherical. Then, 
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once 7 drops below 6 (j l and the jet begins to spread, the flux at 
frequencies above the self absorption would decrease as f" 1//3 . 
At lower frequencies which are in the self absorbed regime the 
flux will be a constant until the self absorption frequency is 
reached. 

These predictions are different from those derived by 
Rhoads, who considered the case when v m < v a where he 
found that the flux rises linearly with time. However, based on 
GRB 970508, it seems that this regime of v m < v a is relevant 
only for very late times, about a hundred days after the burst. 

At high frequencies two light curves are possible, depend- 
ing whether the radiating electrons are cooling {v > v c ) or not 
(v < v c ). The slope itself also depends on the electron power 
low distribution index p. Below the cooling frequency we ob- 
tain 
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Above the cooling frequency we have 
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Note that for a spreading jet, the light curve decay index (but 
not the spectrum) is independent of whether v > v c or v < v c . 
This is due to the fact that v c is constant in time in the case of 
a spreading jet. Since p determines both the light curve and 
spectrum, a parameter free relation between the temporal decay 
index a and the spectral index j3 can be given. For a spherical 
expansion we have: 



a = 
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While for an expanding jet we have: 
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These results are summarized in Table 1 . 

3. OBSERVATIONS 

3.1. GRB 980519 

GRB 980519 was one of the brightest of the bursts detected 
in the BeppoSAX WFC (Muller et al. 1998; in 't Zand et 
al. 1999), second only to the recent GRB 990123 (Feroci et 
al. 1998). The BATSE measured fiuence above 25 keV was 
(2.54 ±0.41) x 10~ 5 ergs cm" 2 , which places it among the top 
12% of BATSE bursts (Connaughton 1998). An X-ray obser- 
vation with the BeppoSAX Narrow Field Instruments detected 
an afterglow (Nicastro et al. 1998). GRB 980519 had the most 
rapid fading of the well-documented GRB afterglows, consis- 
tent with r 2 ()5±0 04 in BVRI (Halpern et al. 1999). The power- 
law decay index of the X-ray afterglow, a x = 2.07 ±0.11 as 
reported by Owens et al. (1998), is consistent with the opti- 
cal. The X-ray temporal decay of GRB 980518 is the fastest 
of the seven afterglows that were well measured by BeppoSAX 
(Owens et al. 1998). The optical spectrum alone is well fitted 
by a power law of the form I /- L20 ± 25 5 while the optical and 



X-ray spectra together are adequately fitted by a single power 

lawi/ -1.05±0.10_ 

The relation between the spectral slope and the temporal de- 
cay is inconsistent with the simple spherical fireball model that 
predicts that the time decay light curve index, a, and the spec- 
tral shape power law index, /3, are related according to equation 



10 This inconsistency is independent of the value of p. These 
observations are consistent with each other if we assume an ex- 
panding jet phase for which equation [Tl| is applicable. 

It is difficult to determine the exact value of p from these 
observations. However, we note that they are consistent with a 
value of p ~ 2.4 that arises in other bursts. This will fit the op- 
tical power law decay if the electrons are not cooling i.e., v c is 
above the optical band. It will also fit the optical spectral index 
which has large uncertainty. The optical to X-rays slope is inter- 
mediate between the value obtained for slow cooling (~ -0.8) 
and that obtained for fast cooling (-1.25). This indicates that 
the cooling frequency is between the optical and x-rays. 

The fact that this transition took place less than 8.5 hours af- 
ter the burst shows that the opening angle of this jet was rather 
small: < 0. 1, leading to a beaming factor of 300 or larger! We 
note that the two strongest GRBs detected by the BeppoSAX 
WFC are inferred to have a large beaming factor. This may in- 
dicate that a significant fraction of the spread in luminosities is 
given by the beaming effect. 

3.2. GRB 990123 

GRB 990123 was a remarkable burst with a very high GRB 
fiuence and with a prompt optical and X-ray emission. We in- 
terpret this emission as resulting from the early forward shock 
(X-ray) and from the early reverse shock (optical). The reverse 
shock has also produced the early radio flare (Sari and Piran 
1999). The forward shock emission is directly related to the 
later optical and X-ray emission, while the reverse shock emis- 
sion decayed like f" 2 , and disappeared quickly. The late opti- 
cal afterglow showed a power law decay with f- 11±0 03 . This 
behavior continued from the first late observation (about 3.5 
hours after the burst) until about 2. 04 ±0.46 days after the burst. 
Then the optical emission began to decline faster (Kulkarni et 
al. 1999). The simplest explanation is that we have observed 
the transition from a spherical like phase to an expanding jet 
phase. The transition took place at ~ 2 days, corresponding to 
6*o ~ 0.1. This implies a beaming factor of about 100 reducing 
the energy of the burst to 3 x 10 52 ergs. This is the first, and by 
now the only, burst in which such a break was detected. The 
decay before the break is well measured and fits an electron 
distribution with p = 2.5. 

3.3. GRB 980326 

GRB 980326 was another burst with a rapid decline. Groot 
et al. (1998) derived a temporal decay slope of a = 2.1 ± 0.13 
and a spectral slope of /3 = 0.66 ± 0.7 in the optical band. Such 
rapid temporal decay suggests a jet like evolution. As Groot et 
al. (1998) note, the large uncertainty in the spectral index al- 
lows in this case also a spherical expansion interpretation (with 
somewhat unusual values p = 4.2 or p = 5.2). However, this 
measured temporal decay was dependent upon a report of a 
host galaxy detection at R = 25.5 ±0.5, which was included 
as a constant term. The detection of a host has since been deter- 
mined to be spurious; better data show no constant component 
to a limiting magnitude of R = 27.3 (Bloom & Kulkarni 1998). 
When the previously assumed constant component is removed, 
the overall light curve is concave, in disagreement with a jet 
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Table 1 

The spectral index p and the light curve index a as function of p. Typical values are quoted using p = 2.4. The 

PARAMETER FREE RELATION BETWEEN a AND j3 IS GIVEN FOR EACH CASE (ELIMINATING p). 



interpretation. If the last detection is interpreted as a different 
phenomenon (Kulkarni, 1999) then the remaining points show 
a rapid decline - in agreement with a jet. 

3.4. GRB 970228 and GRB 970508 

In both GRB970228 and GRB970508 there was no observed 
break in the light curve as long as the afterglow could be ob- 
served. GRB970228 was observed by HST six months later, 
at which point it was still following a power-law decay as 
r i.i4±o.05 ( Fruchter et al. 1998). GRB970508 was observed 
for 9 months to decline as r 1 23±0 04 (Zharikov et al., 1998), at 
which point became as faint as its host galaxy. This set a limit 
on the beaming in these events of 9q > 1 . The beaming factor 
is therefore less than an order of magnitude. 

4. CONCLUSIONS 

We have seen indication of a jet-like behavior in three bursts. 
Two other bursts did not show any break in their optical light 
curves which have been observed for a long time. In several 
other bursts the situation is inconclusive and their short after- 
glow is consistent with rather narrow jets. We suggest that jet- 
like behavior is the common one in GRBs. Moreover, the range 
of possible beaming angles, from 0o < 0.1 for GRB 980519, 
to ~ 0.1 for GRB 990123 to > 1 for GRB 970228 and 
GRB 970508 is quite large. These beaming angles are consis- 
tent with the limits set by searches for "Orphan" radio (Perna 
and Loeb, 1998) and X-ray (Grindlay 1999) afterglows. 

The suggestion that GRBs are beamed has several implica- 
tions. First, this implies that the GRB "inner engines" must 
include a collimation mechanism in addition to the required 



acceleration mechanism. This makes the similarity between 
GRBs and some AGNs, more specifically Blazars, even greater. 
Second the beaming reduces the energy budget of this phe- 
nomenon. Beaming of 0o ~ 0.1 reduces the required energy by 
a factor of 400. Interestingly, the evidence for jets arises most 
clearly in the two strongest bursts detected by BeppoSAX so 
far. It may provide a hint on the energy budget and on the effect 
of beaming on the luminosity function. GRB models based on 
"regular" compact objects become more appealing once more. 

What would be the effects of beaming on the observed lumi- 
nosity function of GRBs and on possible intrinsic correlations 
between different features of the source? One might expect that 
different viewing angles might have a very strong effect on the 
observed luminosities and other characteristics of the GRBs. 
However, as implied by the compactness problem typical ini- 
tial Lorentz factors during the GRB, 70 > 100. An observer 
that is more than 7,7' away from a GRB jet will practically miss 
the GRB. Observers, within this range will see a weaker, softer 
and longer GRB. However, Mao and Yi (1994) point out that if 
0o ^S> 7q 1 then the fraction of bursts that are viewed "sideways" 
is rather small. It would be approximately J^/Oq, a few per- 
cent, if we use the value inferred for GRB 990123. Thus this 
will introduce a small population of intrinsically weaker, softer 
and longer bursts. Since these bursts will be weaker, we expect 
that we would observe only a small fraction of them. 

This research was supported by the US-Israel BSF grant 95- 
328, by the Israeli Space Agency, by NASA grant NAG5-3516 
and by the Sherman Fairchild Foundation. 
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